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^ • We consider the hadronic B decays B'^ — > D^t], B^ -^ D*'^r], B^ — > D^ri\ 

B^ -^ D*^ri' in the framework of a quark-flavour basis and factorization. The 
formahsm allows to compute the decays to the rj meson and to relate them to 
those of the 77'. Measuring the branching ratios of these processes may shed light 
on the nature of the rj-rj' mixing. On the experimental side, only upper limits 
on the branching ratios are known at present. 



I. INTRODUCTION 

The hadronic B decays 5° -^ D^r], 5° -^ D*^t], 5° -^ D^ , 5° -^ D*^i are of 
considerable interest both theoretically and experimentally. Measuring the branching ratios 
of these processes may shed light on the nature of the 77-77' mixing. On the experimental 
side, only limits on these decays are known at present from CLEO |[T| . A theoretical estimate 
of B^ -^ D^*^^r], B^ -^ D^*^^r]' is important also because these processes are a background 
for the semi-inclusive processes B^ —>■ X^t], B^ —>■ Xgi]'. For a recent experimental analysis, 
seei. 

The processes B^ -^ D^*^^r], B^ -^ D^*^^r]' proceed via the internal spectator diagram 

shown in Fig. 1. They are colour suppressed since the colour of the quarks from the virtual 

W has to match the colour of the produced c quark from b —>■ c and the spectator antiquark 

coming from the B meson. Here we consider the quark component rjg of the rj particle in 

the so called quark- flavour basis [§]: 

/r/\ /COS0 -sin0\ /r]g\ 

\r]' J \sin0 COS0 / \r]sJ 

where the Fock state description of r^g ,, is: 

, , ^ \uu + dd) ,„. 

\V,) = ^, ^ ' +■■■ (2) 

\r]s) = ^s|ss) + ... (3) 

The ellipses indicate higher Fock states, like glue states, while the wave functions \E'g,s are 
the light-cone wave functions of the parton states. Our assumption is that the rig component 
is the one to be considered in the factorization diagram resulting from Fig. 1 (shrinking to 
a point the W propagator). We will exploit the relation between 77, rj' and ?7g to compute 
the amplitude for the process B^ -^ D^*^^r]'. 
We recall the form factor parameterization: 



{r,{q,)\V'^{q)\B\p)) 



rni — m 



{p + q,r + ^2 
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with -Fi(O) = -Fo(O). Moreover we need the matrix elements: 



Fi{q' 



Foiq% (4) 



{YAClV^lDip)) = ifop^ (5) 



{VAC\V''\D*{e,p)) = e^niD^fo*, (6) 

and we will use the values [| /d = 200 MeV, fo* = 230 MeV. 

Adopting the Bauer-Stech-Wirbel approach to factorization [^ we can consider the dia- 
gram in Fig. 1 describing the t] coupling through its dd component. The amplitude we need 
to compute is: 

GsoDOr, = (r7D°|i/^^^|5°) = cos0(r7,D°|i7^sw|^o^ 

= ^V.,V:,a,fnFr''\ml){ml - m^f-^, (7) 

where is determined in [|| to be = 39.3° and 02 is the phenomenological coefficient of 
the Bauer-Stech-Wirbel effective Hamiltonian 02 = 0.29. It is significantly larger than the 
leading-order result a\^ ~ 0.12 corresponding to naive factorization. Such a large correction 
is not surprising. The QCD factorization formula of can not be used to compute rigorously 
02 in these decays, however a rough estimate in the limit ttIc -C rrif, gives 02 ~ 0.25 [0]. The 
factor l/-\/2 in (^ accounts for the fact that only the dd component enters in our diagram. 
The strange quark component is absent. 

The amplitude Ggo/jo^/ is then connected to (0) through: 

Grj^o£)0„/ = tan0 Grj^o£)0„. [o) 

In an analogous way we find: 

GsoD*o, = (77D*°|if^^W|5°) = cos0(r7,D*°|i7^sW|^o^ 

= ^V.,V:,a,fn*mn^Fi^-'^\ml,) e ■ {p + qj ^ . (9) 

The term e- [p + qri) is easily evaluated summing over the polarisations oi D* once the square 
modulus of this amplitude is considered. 

II. NS MODEL 

In order to get an idea about the amount of model dependence, we will consider two 
different models for the form factors. We start by considering a model by Neubert and 
Stech (NS) which is a parameterization based on simple assumptions, in reasonable 
agreement with the available data and with most theoretical predictions. In the case of 



heavy-to-light transitions, the form factors no longer obey the symmetry relations valid in 
the heavy-to-heavy decays. One can account for this by introducing, for each form factor, a 
function ^i{w) replacing the Isgur-Wise function: 

Zy/mBmM 
2N 2^mBmM w + 1 
mB + tum 2 

where 

„^.„..,, = '"^ + '"'^-^' . (12) 

ziJiBrnM 

The maximum value of if = Wmax is obtained at q^ = 0. For an estimate of the functions 
C,i{w) a simple pole model is used: 



^iM = \ — — ^ , (13) 

Vw + ll + r w — Wi 



^o(w) = J ^ (14) 

* 1 + r - 

"'max -L 

where 

AmBmv \ M^ - {ms - mvY j 
where V is the vector meson in the same doublet as the scalar meson involved in the decay 

described by the Fq form factor (for the r] this corresponds to uj). Moreover 

ml + ml^ - Ml 
"^ = 2mBmM ' ^'^^ 

where Mi is the mass of the nearest resonance with the appropriate spin-parity quantum 

numbers (for Mq = M^^ = 5.754 GeV it is the 0+ pole and for Mi = Mb^ = 5.325 GeV it 

is the 1~ pole). 

These form factors satisfy the relation 

Fi(0) = Fo(0), (17) 

and scale as ^i{w) ~ w^'^/^ for large w, in accordance with the scaling rules obtained in [Q. 
The B -^ 1] form factors we obtain in this model are given in Table 1. The results for the 
branching ratios are given in Table 2. We also compute the corresponding B ^ n decays 
in order to allow a comparison with experimental data and theoretical predictions that is 
independent of the assumptions concerning the tj mixing. 



III. CQM MODEL 

The computation of the form factor Fq and Fi can be carried out with the aid of a Con- 
stituent Quark Model (CQM) [^. This model has shown to be particularly suitable for the 
study of heavy meson decays. Since its Lagrangian describes Feynman rules for the ver- 
tices (heavy meson)-(heavy quark)-(light quark) [|I^, transition amplitudes are computable 
via simple constituent quark loop diagrams in which mesons appear as external legs. We 
determine Fq ^^^ considering the 1]^ as a Goldstone boson. Therefore one has to take into 
account the relations between rjg and rjq for the process under study. To begin we observe 
that: 

q^iVslV^m = (m| - m^)Fi''-^'\ (18) 

where q'^ is the momentum carried by the current V^. For our purposes rj is essentially the 
Goldstone octet component rjs'- 

{r]\V^\B°) ~ {r]s\V^'\B^). (19) 

On the other hand the amplitude related to the process in Fig.l is: 

{r]\V^'\B'^) =cos(j){r]g\V''\B^), (20) 

if we also take into account a small mixing angle, 6s in the notations of [H], we write: 



{r^slV^lB') = ^{v,\V^m. (21) 



Again we stress that here we are neglecting the singlet component and we consider that the 
B ^ 1] matrix element which enters in the factorization diagram is the one containing only 
the dd components. In other terms we can write: 



q'{Vs\V^\B') = ^« - K)^^ , (22) 

in such a way that: 

p{B^,,) _ cos Os (m| - ml) .^^^^^ 
COS0 (m^ — mf^^) 

Notice that ds is determined in to be 6s = —21.0°. For the Fi form factor we obtain: 

^' ^^^*)-6-(p + gJcos0^^ ^"'^*^- (^^^ 
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We compute the tj particle via its r^g component with a Goldstone hke couphng in the 
meson-quark loop computation. The interaction Lagrangian generating this coupling has 
the form: 

i'al ■ Aabl^i'b, (25) 

where the indices a,b are the u,d,s quark indices, ip being a triplet of flavour- S'f/3. The 
structure of the A matrix is: 



A, = -^m^^-^%0, 



(26) 



and ^ = e*'^/-^. The n matrix has the well known form: 



V2^ VE 



n 



\ 



K- 



1T 
'V2~ 



m 

V6 



K+ 



K^ -\ 



(27) 



TS^s/ 



When the dd component is taken into account the Feynman rule for the vertex is 
Q^lfilb/iV^fTr)- We isolate three contributions to the Fqi form factors: a non-derivative 



contribution, a direct and a polar contribution |11]; see Fig. 2. The direct and the non- 
derivative contributions are represented in the CQM model by loop diagrams in which the 
internal lines are heavy quark and light quark propagators, while the external legs contain 
the incoming meson and the vector current. In the non-derivative diagram the rj external 
leg is attached at the same vertex (heavy quark)-(light quark)-(current) due to the structure 
of the CQM interaction Lagrangian. In the direct diagram the 77 is instead attached to the 
light quark internal line (see fill). The polar diagram allows for an intermediate polar state 
separating the current insertion from the loop diagram. Each of these diagrams contributes 
to the form factors Fqi in a calculable way. Summing up these contributions and imposing 
the condition -Fi(O) = -Fo(O), which eliminates the spurious singularity in the form factor 
decomposition, one obtains: 

F,{q') = aM'){Ff + Ft{q') + Fr\q')), (28) 

where z = 0, 1. The functions aj(g^), j = 0, 1 are needed to impose the condition -Fi(O) = 
Fo(0). This is because the sum in parenthesis reproduces the form factor on the Ihs of the 
previous equation only at g^^^, i.e., a(gLx) 



1. Following |]Tl| we choose: 



aj{q^ 



1 + a,- 



m 



B 



mz 



2mBAy 



1 - ajErJK 



X' 



(29) 



where A^ = 1 GeV. In such a way we enforce that a(q'max) — 1- ^r^ ^^ ^^^ V energy in the B 



rest frame. We borrow oq from [O where we had Oq = 0.27. The form factor used in the 



computation of the branching ratios are of course those given in (^). The non-derivative 
(nd) contributions to the form factors are given by: 



nnd JB T:^nd JB 



pnd _ J^ pnd ^ _i^. (30) 



The direct contributions are: 



^Dir(^2) ^ S^ + S^ ^^^^ 



2 

2A{q^ 



F^il') = ^^-^' + 'lf\^^^ + ^'\ (32) 



where: 



and: 



. ^ (J<^^ (33) 



.2 /<-)„„ ^2 
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In these expressions: 



" ■ ^^ = — t;^^ — • (^^) 



The symbols Z, Zj^, fij indicate integral expression listed in [|1^. These are functions of A//, 
the difference between the heavy meson and the heavy quark constituent mass mu — m-g, 
A = A// — V ■ Qri and v ■ qr/, v being the four velocity of the B meson, g^ is the momentum 
carried by the current insertion in the CQM loop diagram, m is the constituent light quark 
mass, m = 300 MeV. The suffix H is referred to the H = (0^, 1^) multiplet predicted by 
heavy quark effective theory. 
The polar contributions are: 



,poi. 2n ^ 1 (hrr^^/m^F+\ 



FrW = -Z^--^. 'v.. . .2.^2 . (36) 



m 



B-^vK v^/^ / 1 " 1^/^ 
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where B** is the 0+ state of the multiplet S = (0+, 1+) while B* is the 1~ state of the H 
multiplet. h is the couphng constant HSir and in CQM one finds h = —0.76 ± 0.13 while 
g = 0.46 ± 0.04, from HHtt. We use the central values for all these constants and the value 
Ah = 0.4 GeV. Variations of A^ in the range 0.3 — 0.5 GeV induce few % variations in the 
form factors (see Table 2). F and F^ are defined as the decay constants of the H and 5* 
heavy mesons respectively. 

Note that the form factors calculated in this section already include the factor taking into 
account that the process we are considering only involves the dd component of t] through 
the numerical coefficient of the Feynman rule. Therefore when combining (0) and p3| ) one 
has to be careful to remove a factor l/-\/2 from (|^): 

in terms of the t]s state. From (^) and (P) we get in a similar way: 



GbodOt, = ^VcbV*aa2fDFi^^'^'^ {ml){ml - m^) cos Og (37) 



GsoD'Ori = -^VcbV*^a2fD*mD-F[^~''^''\m\),)e ■ {p + q^) cos6'8 . (38) 

Numerical results are given in Table 1 and 2. The two models give similar results for the 
form factors and branching ratios. The amount of model dependence due to the form factors 
is difficult to quantify. By comparing the two models one can conclude that they affect the 
branching ratios at the level of ~ ±0.1 x 10~^. The a2 coefficient is considered fixed at 
the value 0.29 in the calculation. As it enters as a^ in the branching ratio, by changing its 
value to 0.25 may add to the error on the branching ratio up to ~ ±0.2 x 10^^. However if 
the quark picture used in the calculation is correct such a large variation is not expected as 
B'^ -^ D^i] is in this case on the same footing as B'^ -^ D^tt. 

IV. CONCLUSION 

In order to test the idea of relating the rj and rj' using the quark-flavour basis in the decays 
considered in this work, experimental data will have to measure the ratio of Br's i]'/ri. If 
different from what estimated, it probably means that the quark mechanism (see Fig. 1) 
is not sufficient to explain the decays into rj'. This would suggest to look at some other 
mechanism, like the gluon anomaly explored in [O] where the large observed production 
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of 7]' m. B ^ Xgi]' decays has been studied. The mechanism suggested there is based on 
the subprocess b — > sg* -^ srj'g where the virtual gluon emerging from the standard model 
penguin couples to rj' via a gluon anomaly vertex g*gr]'. The possibility that the other gluon 
g is emitted by the light quark inside a B meson has been examined in W^. The B -^ Krj' 



decay has been recently examined in the context of perturbative QCD in |1J|. In the case 
oi B ^ Drj' studied here it is more difficult to imagine some gluonic mechanism of the 
kind described before. Anyway we cannot exclude this possibility since the t]' can have a 
large glue component [J13l; see also [|15]. Under the assumptions considered in this work 



B^ -^ D^f] is similar to B^ -^ D^tx. The calculation is related to the one for B^ -^ D^rj' 
using the quark-flavour basis. We have checked that two different ways of computing the 
form factors give similar results. 

Note added in proofs 



Very recently there was a phenomenological determination of the 02 coefficient |16| based 
on new data from Belle |T^. This seems to be in contrast with all previous determinations 
based on the factorization approach addressing the problem of final-state interaction in 
the B -^ Dn channel. Such a problem could arise also in this analysis where the B -^ Dr] 
is considered as the starting point, and the factorization value for 02 is assumed. We will 
consider this aspect in future work. 
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TABLE I: Values of the form factors used in the calculation. These are used in slightly different 
formulas for computing the amplitudes. The CQM form factors are indeed extracted using a 
Goldstone-like coupling for the rj particle and refer only to the dd component, while those of the 
NS model refer to the uu + dd (see Sec. 3). In order to compare the two models on the same footing 
the form factors of the CQM model have to be multiplied by v2. 



F^-\Q) F^-\ral) F^-\ral.) 



CQM Model 0.15 0.17 0.23 

NS Model 0.26 0.28 0.33 



TABLE II: Theoretical predictions and 90% C.L. upper limits on branching ratios. The errors 
quoted for the CQM model only refer to the variation of the parameters in the model and not to 
the model dependence of the result. 



Decay Mode 



C.L. m NS Model CQM Model 



B^ 



D^TT < 1.2x10-^ 0.77 X 10"^ 1.3+^1 x lO""^ 



0.3 



D*^TT < 4.4x10"^ 1.05 X 10"^ 1.1 ±0.3 X 10"'' 



B° ^ D^T] < 1.3x10"^ 0.50 X 10"^ 0.44 ±0.02 X 10"*^ 

^0 ^ 2)*°77 < 2.6x 10-^ 0.60 X 10"^ 0.70 ± 0.04 x 10"*^ 

B° -^ D^r]' < 9.4x 10-^ 0.32 x 10"^ 0.30 ± 0.02 x 10"^ 

B° -^ D*°r]' < 14x10"'' 0.41 x 10"^ 0.47 ± 0.04 x 10"^ 



D< 




u 



BO 



W" 



d 



Tic 



FIG. 1: Flavour diagram for B^ -^ viv')- The r/g component allows to bridge between the two processes. 
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FIG. 2: CQM diagrams. The heavy line represents the incoming heavy meson. The double line is the 
heavy quark, the curly line represents the current insertion, (a) is the non- derivative diagram in which the 
T] is coupled at the same vertex as the current, (b) is the direct diagram and (c) is the polar diagram: an 
intermediate heavy meson state is taken into account. 



